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Biological cells can be treated as composites of graded material inclusions. In addition to biomaterials, graded 
composites are important in more traditional materials science. In this article, we investigate the electrorotation 
(ER) spectrum of a graded colloidal suspension in an attempt to discuss its dielectric properties. For that, we 
use the recently obtained differential effective dipole approximation (DEDA) and generalize it for non-spherical 
particles. We find that variations in the conductivity profile may make the characteristic frequency red-shifted 
and have also an effect on the rotation peak. On the other hand, variations in the dielectric profile may enhance 
the rotation peak, but do not have any significant effect on the characteristic frequency. In the end, we apply our 
theory to fit experimental data obtained for yeast cells and find good agreement. 

PACS numbers: 82.70.-y,77.22.Gm,77.22.-d,77.84.Lf 
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I. INTRODUCTION 

AC electrokinetic phenomena, i.e., electrorotation, dielec- 
trophoresis and traveling wave dielectrophoresis, have gained 
an increasing amount of attention during the past decade. This 
is due to their wide range of applications from cancer research 
to identifying and separating parasites, cell populations and 
viruses [[y, 0], and even to design of nanomotors [|| H]. De- 
spite the number of applications, there is a need for a theory 
that treats the different aspects of electrokinetic phenomena on 
an equal footing starting from the general underlying physical 
principles 

In this article we concentrate on electrorotation. Electro- 
rotation (ER) is a phenomenon in which an interaction be- 
tween a rotating ac electric field [g |7|, || and suspended 
dielectric particles leads to a rotational motion of the par- 
ticles. This phenomenon was first observed experimentally 
over four decades ago [|| and during the last two decades, ER 
has been increasingly employed as a sensitive tool for non- 
invasive studies of a broad variety of microparticles, ranging 
from living cells to spores, seeds as well as synthetic materi- 
als[|00||||0|. 

In electrorotation, a rotating ac electric field induces a 
dipole moment inside a polarizable particle which rotates at 
the angular frequency of the external field. If the frequency of 
the field increases, the period of the rotating field becomes 
comparable to the time scale related to the formation of a 
dipole. To minimize its energy, the dipole tries to line up with 
the field but at high frequencies the particle is no longer able to 



follow the field. This leads to a lag in its response and a torque 
is induced by the interaction of the out-of -phase part of the in- 
duced dipole moment and the external field. The torque and, 
in turn, the rotating speed are proportional to the imaginary 
part of the dipole factor (also called Clausius-Mossotti factor). 
In addition, the rotating speed of a particle is inversely propor- 
tional to the dynamic viscosity, since the particle is suspended 
in a viscous medium and experiences a dissipative drag force. 

Different models have been suggested and used for analyz- 
ing experimental results, see e.g. Refs. [g [B|] and refer ences 
therein. The most commonly used ones are the so-called shell 
models. This approach was first used by Fricke already in 
1925 jl6[]. A model for a biological cell consisting of a single 
spherical dielectric shell was introduced by Arnold and Zim- 
mermann in 1982 to investigate the rotation of a mesophyll 
protoplast [^]. Four years later, a two-layer model consisting 
of two spherical shells is put forth to discuss the electrorota- 
tion of a single plant cell jl7|]. Due to the inhomogeneous 
compartmentalization of cells, three-shell models have been 
used to investigate ER of liposomes [B]. 

These inhomogeneities can be studied in the framework of 
graded materials, i.e., materials with spatial gradients in their 
structure. Biomaterials such as cells and bones are typical ex- 
amples as their composition varies through the object. Grada- 
tions can also be used to control and improve the strength and 
other properties, for a recent review see Ref. [ 18 1. In biolog- 



ical and medical applications identification and separation of 
graded objects in micron scale (and below) is of great impor- 
tance. In these situations gradation is generic due to inhomo- 
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geneous compartments of cells, the cytoplasm often contains 
organelles or vacuoles. These inhomogeneities are reflected 
in their dielectric properties. Since differences and changes in 
the dielectric properties are intimately related to physical and 
chemical properties of cells (as well as other objects), electro- 
rotation and other electrokinetic phenomena allow identifica- 
tion of particles and detection of compositional changes. 

As far as one decade ago, Freyria et al. observed a graded 
cell response when they experimentally study cell-implant in- 
teractions [|l9|]. Recently, motivated by positional informa- 
tion of protein molecules, Honda and Mochizuki performed 
computer simulations of cell pattern formation showing that 
graded cells can be generated by cell behavior involving 
differences in intercellular repulsion p0|]. To discuss such 
graded cells, the above-mentioned multi-shell models fail to 
apply. It is thus necessary to develop a new theory to study 
the effective properties of graded composite materials under 
externally applied fields. In a recent work we succeeded 
in deriving the dipole factor after putting forth a differential 
effective dipole approximation (DEDA). In this work, we gen- 
eralize DEDA to non-spherical particles and apply it to inves- 
tigate the ER spectrum of graded particles in an attempt to 
investigate their dielectric behavior. Finally, to demonstrate 
the validity of our model, we fit our theory to experimental 
data [22 1 and find good agreement. 



II. FORMALISM 

We consider an inhomogeneous biological cell or a parti- 
cle of radius a with complex dielectric constant profile 1\ (r), 
embedded in a host fluid of dielectric constant e 2 . Here 
e = e + <r/(2iTif), where e denotes real part of the dielec- 
tric constant, a stands for conductivity, / is the frequency of 
an external field, and i = V~ T- The dipole factor reflects the 
polarization of a particle in surrounding medium. 

In a recent work J2l| ] , we derived the dipole factor for spher- 
ical particles by introducing a differential effective dipole ap- 
proximation (DEDA). The idea of DEDA can be summarized 
as follows: Consider a shell model for an inhomogeneous par- 
ticle. In DEDA one adds new shells of infinitesimal thickness 
to the particle under consideration. Each of these cells have 
distance dependent dielectric constant. As the thickness of 
the layer approaches zero, i.e., dr — > 0, then the correction to 
the dipole factor due to it is infinitesimal and one can obtain 
a differential equation for it - hence the name differential ef- 
fective dipole approximation. Next we extend DEDA to the 
general case of graded spheroidal particles. Then, we study 
the effect of different conductivity and dielectric profiles to 
the electrorotation spectrum, and finally compare the theory 
to experiments. 

We first consider a homogeneous dielectric spheroid in a 
uniform electric field. The spheroidal inclusion has a dielec- 
tric constant ii and it is embedded in a host medium of di- 
electric constant e 2 . A spheroidal particle can be parametrized 



by x 2 + y 2 



{z/hf 



and the dielectric profile for the 



graded particle is given by e(r), where < r < a' and h 
is the aspect ratio and a' is the length of the semi principle 



axis along x — (y) axis. It is worth noting that h > 1 de- 
notes a prolate spheroid, whereas h < 1 indicates an oblate 
one. We assume that the longest axis of the spheroid is along 
the z— axis with the depolarization factor L z , satisfying a sum 
rule L z + 2L xy = 1 where L x r y ) is the depolarization fac- 
tor along x(y)~ axis. For a homogeneous spheroid, the dipole 
factor along z— axis b z is given by 



ei - £2 



3 h + L z {h - hY 



(1) 



which measures the degree of polarization for a particle in 
an external field. Next, we add a confocal spheroidal shell 
of dielectric constant e to the spheroidal particle to make a 
coated spheroid. Then, the dipole factor of the coated spheroid 
is[|3}0] 



h z 



1 (e-e 2 )+3a:ip[e + Ls(e2-e)] 

3 [h + L z (e - h)] + 3xi P L z (l - L z )(e - e 2 ) 



(2) 



where p is the volume ratio of the core with respect to the 
whole spheroid, and 



ei 



3 e + L z (ii - e) 



This approach can be extended to more shells of different 
dielectric constants, at the expense of obtaining more compli- 
cated expressions. It is easy to check that b\ z (Eq. (Q)) reduces 
to b z (Eq. ([[])) when e = ej.. Thus, the dipole factor remains 
unchanged if one adds a spheroidal shell of the same dielectric 
constant. 

Next, we consider an inhomogeneous spheroid with the di- 
electric profile e(r). To establish DEDA, we mimic the graded 
profile by a multi-shell construction by building up a dielectric 
profile by gradually adding shells. We start with an infinitesi- 
mal spheroidal core of dielectric constant e(0) and continue to 
add confocal spheroidal shells of dielectric constant e(r) until 
r = a' is reached. At r, we have an inhomogeneous spheroid. 
To establish a functional form, we replace the inhomogeneous 
spheroid by a homogeneous spheroid of the same dipole fac- 
tor, and the graded profile is replaced by an effective dielectric 
constant e(r). Thus, Eq. (jj]) becomes 



C-2 



3 e 2 + £,(e(r) - e a ) 



(3) 



Next, we add a confocal spheroidal shell of infinitesimal 
thickness dr, and dielectric constant e(r). The dipole factor 
will change according to Eq.(^|). The effective dielectric con- 
stant e(r), being related to b z {r), will also change. To see this, 
let us write b\ z = b z + db z and take the limit dr — > 0. We 
obtain a differential equation: 



dr 



[e 2 (l + 3b z L lz ) - (1 - 3L z b z )e(r)} 



re 2 e(r) 

[e 2 L z (l + 3b z L lz ) + L lz (l - 3b z L z )~e(r)] , (4) 



where L lz — 1—L Z , The dipole factor of a graded spheroidal 
particle can be calculated by solving the above differential 



3 



equation with a given gradation profile e(r). This nonlinear 
first-order differential equation can be integrated, at least nu- 
merically, if we are given a profile for e(r) and an initial con- 
dition for b z (r = 0). It is worth noting that the dipole factor 
is a tensorial quantity and the dipole factor along x(y)— axis 
b x ( y ) can be obtained by changing L z to L x ^ in the equation 
above. It is important to note that substituting L z = 1/3 into 
this equation yields the DEDA equation for spherical graded 
particles [^TJ, 



volume average: 



db 
dr 



[(l + 26)e 2 -(l-6)e 1 (r)] (5) 



3re 2 ei(r) 
[(l + 26)e a + 2(l-&)ei(r)] 



where < r < a, and a is the radius of the spherical particle. 
This is the expression we will use in our studies of the effect 
of conductivity and dielectric profiles on the ER spectrum. 

The dipole factor of a graded spherical particle can be cal- 
culated by solving the above nonlinear first-order differential 
equation (Eq. (Q)) with a given profile ?i(r) and the initial 
condition 



b(r = 0) 



ei(0)-g 2 
e!(0) + 2e 2 



(6) 



If we apply a rotating field to the system, a torque be- 
tween induced dipole moment inside the cell and the elec- 
tric field vector will set the cell in a rotational motion. It 
is well known that the frequency-dependent rotation speed 
is proportional to the imaginary part of the dipole factor 
— Im[6(r = a)]. The rotation speed is given by 



«(/) 



•— — -Im[&(r 
2r] 



(7) 



where E Q is the magnitude of the rotating field and rj the vis- 
cosity of the host fluid. Hence, when b(r — a) is solved 
by integrating Eq. (||), we can investigate the electrorotation 
spectrum. 

In order to perform numerical calculations, we take the con- 
ductivity and dielectric profiles to be 

a x (r) = CTl (0)(r/a) n , r<a 
ei(r) = ei(0) + c(r/a), r<a 

where n and c are profile dependent constants. The profile is 
quite physical in the sense that the conductivity can change 
rapidly near the boundary of cell and a power-law profile pre- 
vails. The profile constant n can take any positive values and 
it can be larger than unity. From the above equations it is 
clear that the conductivity is largest at the boundary. This 
is the case in real systems as well, since the outermost shell 
(e.g. cell wall) of a cell is always conductive and often has 
the largest conductivity. On the other hand, the dielectric con- 
stant may vary only slightly and thus a linear profile suffices. 
In particular, the dielectric constant at the center, i.e., ei(0), 
may be larger than that at the boundary. Thus, we will choose 
c < 0. By integrating the dielectric profile, we obtain an aver- 
age dielectric constant e av for different values of c by using a 



f a 6i(r)r 2 dr 
Jo r2dr 



(8) 



For the present dielectric profile, we obtain ei av — £i(0) 
3c/4. 



III. NUMERICAL RESULTS 

We are now in a position to perform numerical calculations. 
We set o-i(O) = 2.8 x 1(T 2 S/m, ei(0) = 75e , e 2 = 80e , 
er 2 = 2.8 x 1CP 4 S/m for all of the computations. Parameter 
eo denotes the dielectric constant of vacuum. Here, the numer- 
ical results and the comparison to experimental data is limited 
to spherical cells only since we were not able to find appro- 
priate experimental data describing non-spherical cells. For 
the numerical computations, we used the fourth-order Runge- 
Kutta algorithm with step size h — 0.01 to solve the above 
nonlinear first-order differential equation (Eq. (^])). It was ver- 
ified that this step size guarantees accurate numerics. 
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FIG. 1: ER spectrum for profile constants n at c = and c = — 30eo, 
respectively. The spectrum is given as the imaginary part of the 
dipole factor. 

In Fig. U, the ER spectrum is investigated for various in- 
dices n with profile parameters c = and c = — 30eo- 
This corresponds to t\ av = 75eo and 52.5eo calculated using 
Eq. (^]), respectively. Figure [j] shows that the rotation peak 
is always present. It is important to note that the character- 
istic frequency at which the rotation speed reaches its maxi- 
mum may be red-shifted, especially at large n (e.g. n = 5) 
in comparison with small n. Thus, the conclusion is that large 
fluctuations in conductivity of graded particles may lead to a 
red-shift of the characteristic frequency. 
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FIG. 2: Same as in Fig.l, but for various c at n — and n = 1, 



respectively. 



In Fig. g, ER spectrum is investigated for n = and n = 
1, using various slopes c. We use c = 0, c = — 30eo and 
c = — 60eo corresponding to ei av — 75cq, 52.5eo and 30eo, 
respectively. From the figure, it is evident that large spatial 
fluctuations in the dielectric constant may enhance the rotation 
peak. 

From Figs. [I] and |, we find that spatial conductivity fluctu- 
ations are able to make the characteristic frequency red-shifted 
and have a clear effect on the rotation peak. In contrast, spatial 
dielectric constant fluctuations may enhance the peak value, 
but their effect on the characteristic frequency is very small. 
To summarize, we can theoretically obtain a typical ER spec- 
trum with lower characteristic frequency and high rotation 
speed by adjusting conductivity and dielectric profiles. This 
corresponds to the physical/experimental situation in the case 
of graded particles. 

In Figs. [i]and |[ the particles under consideration are spher- 
ical, i.e., L z — 1/3. We should remark here that for non- 
spherical particles, e.g. prolate < L z < 1/3 and oblate 
spheroids 1/3 < L z < 1, the same conclusions hold (no fig- 
ures shown here). In addition, it is also found that for prolate 
(oblate) graded particles, smaller (larger) depolarization fac- 
tor (L z ) can make the characteristic frequency significantly 
red-shifted. 

Finally, to demonstrate our theory, we fit our theory to ex- 
perimental data in Fig. [| The data for yeast cells is extracted 
from the experiments of Reichle et al. [^2|]. In Fig. |J from 
upper panel to lower panel, the symbols denote experimental 
data of the ER spectrum of a yeast cell in water measured at 
the beginning of the experiment, after 20 minutes, and after 
40 minutes, respectively. The corresponding continuous lines 
are the fitted result using the present theory. As discussed 
earlier, the assumption of homogeneous compartments is not 
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FIG. 3: Fitting of experimental ER spectrum of a yeast cell. The 
data is extracted from Ref. p^]. The symbols denote experimental 
data which is measured in water at the beginning, after 20 min, and 
after 40 min, respectively. The corresponding line denotes the results 
using the present model. 



generally justified for biological cells, and the cytoplasm of- 
ten contains organelles or vacuoles. For this reason, the exper- 
imental data is typically fitted using a (two-)shell model. To 
improve fitting using cell models, one would have to adjust 
six parameters, i.e., real dielectric constant and conductivity 
of two shells, as well as cell interior. In contrast, using on our 
theory, only two parameters need to be adjusted, namely the 
profile parameters c and n. 

For all of the three fittings, we set 62 = 80eo, 02 = 



5 



1.1 x 1CT 3 S/m (from the experiment), ei(0) = 75e , cti(O) = 
0.28 S/m, E = 10 kV/m (from experiments), 77 = 2.082 x 
10~ 2 kg/(ms), and a — 5/xm (given by experiment). For 
the upper panel, c = — 30eo (namely t\ av = 52.5eo) and 
n = 0.5. For the middle panel, c = — 20eo (ei av — 60eo) and 
n = 0.8. For the lower panel, c = -25e (ex av = 56.25e ) 
and n — 0.05. 

As seen in Fig. |], our theory performs very well. Due to the 
inhomogeneous compartments of a yeast cell, we believe that 
the graded cell model is a very choice. 



IV. DISCUSSION AND CONCLUSIONS 

In this article, we have generalized the differential effec- 
tive dipole approximation (DEDA) to non-spherical particles. 
This allows studies of shape effects, and indeed colloidal sus- 
pensions and many kinds of cells have non-spherical shapes, 
e.g. red blood cells are oblate spheroids. In particular, for 
oblate spheroids larger depolarization factor leads to a red- 
shift of the characteristic frequency, for prolate spheroids the 
systems behaves the opposite way [14|. To check the valid- 
ity of DEDA, it was compared against exact solutions which 
can be obtained for a power-law profile and a linear pro- 
file [ p6| by solving the Laplace equation for the local electric 
field. In both cases DEDA produced the same results thus 
showing that it is essentially exact. 

In this article, we have applied DEDA to discuss the elec- 
trorotation spectrum of graded colloidal suspension in an at- 



tempt to discuss its dielectric behavior. As a result, we find 
that changes in the conductivity profile can affect both the 
characteristic frequency (red-shift) and rotation peak, whereas 
changes in the dielectric profile may enhance the rotation peak 
only. Consequently, we may conclude that the graded behav- 
ior does play an important role in the dielectric properties of 
graded suspensions. We used our theory (for spherical cells) 
to fit results obtained by Reichle et al. [ ]22| ] for yeast cells and 
obtained excellent agreement. It would be interesting to com- 
pare this theory to experimental results for non-spherical cells 
as well, and hopefully new experiments will be performed in 
that direction. 

Our work may be extended to discuss the ER spectrum of a 
pair of particles by taking into account the multiple image ef- 
fect [ 27 28|] and to investigate high concentration suspensions 
by discussing the local-field effects [29, 3(J. A generalization 
which includes both of these effects is of particular interest. 
In addition to ER spectrum, our theory may be used to discuss 
other electrokinetic effects, e.g dielectrophoresis [ |28[ ^T| ] and 
electro-orientation [32|. 
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